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We study resonant tunneling through a superconducting double barrier structure in graphene 
as a function of the system parameters. At each barrier, due to the proximity effect, an incident 
electron can either reflect as an electron or a hole (specular as well as retro Andreev reflection in 
graphene). Similarly, transport across the barriers can occur via electrons as well as via the crossed 
(specular and/or retro) Andreev channel, where a hole is transmitted nonlocally to the other lead. 
In this geometry, in the subgap regime, we find resonant suppression of Andreev reflection at certain 
energies, due to the formation of Andreev bound levels between the two superconducting barriers, 
where the transmission probability T for electrons incident on the double barrier structure becomes 
unity. The evolution of the transport through the superconducting double barrier geometry as a 
function of the incident energy for various angles of incidence shows the damping of the resonance 
as normal reflection between the barriers increases. 



PACS numbers: 73.23.-b,72.80.Vp,74.45.+c 
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I. INTRODUCTION 

The discovery of graphene, a two dimensional single 
layer of graphite, by K. S. Novoselov et air a few years 
ago, has led to an upsurge in the study of its transport 
properties, both theoretically and experimentally^ "*. In 
graphene, there are six discrete points at the edges of the 
hexagonal Brillouin zone where the energy bands touch 
the Fermi energy, out of which, only two are inequiv- 
alent and are commonly known as the K and K' val- 
leys. The low energy quasiparticle excitations about the 
K and K' valleys behave like massless relativistic Dirac 
fermions. The presence of such quasiparticles in graphene 
provide us with an experimental test bed for observing 
many well-known phenomena in relativistic quantum me- 
chanics, such as the Klein paradojP at low energies. 

The existence of Dirac-like quasiparticles in graphene 
has also motivated a lot of research work in exploring 
the effects due to the proximity of a superconductor. 
Graphene is not a natural superconductor by itself. How- 
ever superconductivity in a graphene layer can be in- 
duced in the presence of a superconducting electrode near 
it via the proximity effeclP. A direct manifestation of 
proximity effect is the phenomenon of Andreev reflec- 
tion (AR)^ in which an electron like quasi-particle inci- 
dent on a normal metal— superconductor (NS) interface 
is reflected back as a hole along with the transfer of two 
electrons into the superconductor as a Cooper pair. Re- 
cently it has been predicted that a graphene NS junction, 
due to the presence of the Dirac-like energy spectrum of 
its quasiparticles, can exhibit specular Andreev reflec- 
tion (SAR)^ ^, in addition to the usual retroreflection ob- 
served in conventional NS junctionsiSl. The presence of 
the SAR process in graphene leads to qualitativ ely dif- 
ferent behaviour i n the tunneling conductance^SQll^ the 
Josephson currenll^5Il51 and the spin current as com- 



pared to those in conventional superconducting hybrid 
junctions. 

An even more intriguing example where the proxim- 
ity effect manifests itself is the phenomenon of crossed 
Andreev reflection (CAR) which can only take place in 
a normal metal— superconductor— normal metal (NSN) 
junction, provided the distance between the two nor- 
mal metals is less than or equal to the phase coherence 
length of the superconductor. This is a nonlocal pro- 
cess where an incident electron from one of the normal 
leads pairs up with an electron from the other lead to 
form a Cooper pair and jumps into the superconductor. 
Due to the presence of the Dirac-like energy spectrum, 
like SAR , graphene can also exhibit specular crossed An- 
dreev reflection (SCAR) in a proximity induced graphene 
NSN junction. The effect of CAR in graphene has been 
studied earlier in Refs. I18I19I in the context of detecting 
entangled states in graphene. However, transport proper- 
ties of a superconducting double barrier (SDB) geometry 
in graphene, i.e. graphene NSNSN junctions, have not 
been studied so far, where resonance effects can be more 
important. 

In general, electronic confinement in graphene is exper- 
imentally challenging due to the effect of Klein tunneling. 
In the recent past, resonant tunneling has been consid- 
ered in doped graphene (single barrier) n-p-n j unctions^ 
and in normal double barrier structurea^IES where res- 
onance effects on the transmission have been investi- 
gated. But here again, the problem of resonant trans- 
mission through doped double barriers in graphene has 
not been investigated. Motivated by this, in this article, 
we study resonant tunneling through a SDB structure in 
graphene, which, because of the correspondence between 
AR and Klein tunneling 1221 would also be valid for a doped 
graphene n-p-n-p-n junction. 

The paper is organized as follows. In Sec. [llj we de- 
scribe the set-up for our system, where two supercon- 




FIG. 1: (Color online) Cartoon of the SDB structure in 
a graphene sheet. Two patches at the two places on the 
graphene sheet depict superconducting material deposited on 
top of it. The schematic of the potential profile seen by an 
incident electron is shown below. 



ducting patches have been deposited on top of a clean 
graphene sheet to form the graphene SDB structure and 
discuss its theoretical modeling. In Sec. |III[ we obtain 
the resonance condition analytically (the Andreev bound 
state condition) by considering multiple (retro) Andreev 
reflections between the barriers. Then in Sec. |IV[ we dis- 
cuss the effect of normal reflection and SAR between the 
barriers on the resonances for angles of incidence other 
than normal incidence. In Sec. |VJ we discuss the numeri- 
cal results of our study and show how the resonances be- 
have as a function of the various parameters in the theory 
and also how the resonances evolve as a function of the 
energy for various values of the incident angle. Finally 
in Sec. |VI[ we present our summary and conclusions. 



GRAPHENE SDB STRUCTURE 



II. 



In our analysis, we consider a clean graphene sheet 
occupying the x-y plane. The SDB structure is formed 
by depositing thin strips of superconducting material on 
top of the graphene sheet at two places. This induces a 
finite superconducting gap (Aie*"^') in the barrier regions 
as a result of the proximity effect of the superconducting 
patches. Here A^ and (pi are the pair potentials and order 
parameter phases on the two patches respectively (i refers 
to the index of the strips). The geometry is shown in 
Fig. [l] The space dependence of the order parameter 
(which also acts as a scattering potential for the incident 
electron) can be expressed as 

V{x) = Ae''^e(a;)e(-a; + a) + Ae"^ 

Q[x - {a + L)]Q[-x + {2a + L)] (1) 

where a is the width of the superconducting barrier in 
graphene and L is the distance between the two barriers. 



Here we assume that the spatial variation of potential 
steps is slow on the scale of the lattice spacing so that 
inter- valley scattering is suppressed. Also here is the 
Heaviside G-function, and we have taken 4>i = (j)2 = </>, 
since we will not be looking at supercurrents (Josephson 
effect) in this work. 

If the width a of the superconducting strips is of the 
order of the phase coherence length of the superconduc- 
tors, the normally incident electron can be transmitted 
across the barriers both as an electron (electron cotun- 
neling (CT)) and as a hole, via the retro crossed An- 
dreev reflection (CAR) process which is nonlocal and oc- 
curs between the same bands. Importantly, besides local 
retro AR and nonlocal retro CAR, graphene can also ex- 
hibit local specular Andreev reflection (SAR) and nonlo- 
cal specular crossed Andreev reflection (SCAR)^^ which 
can switch the valley bands. In the process of retro CAR, 
an electron incident on the barriers is transmitted as a 
hole which retraces the trajectory of the incident elec- 
tron. On the other hand, in SCAR, a hole is transmitted 
nonlocally but in a specular fashion. Here, we restrict 
ourselves to spin singlet (s-wave) superconductors so that 
the electron and the hole are taken from opposite spin 
bands in order to allow the Cooper pair to jump into the 
superconductor with net spin zero. 

Because of the pseudo spin and valley degeneracies 
present in graphene, it suffices to use a four dimen- 
sional version of the Dirac-Boguliobov-de Gennes equa- 
tion (DBDG)^ for electrons and holes which is given by 
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where, U = U(r)+Ep, and the energy e is measured from 
the Fermi level of the superconductor. We assume that 
U{r) = in the normal graphene region and U{r) — Uq, 
a constant, independent of r in the proximity induced 
superconducting region. Note that we have defined di- 
mensionless variables 
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to replace the original ones. 

The solution of the DBDG equations^, describing elec- 
trons and holes with incident energy e inside the normal 
graphene regions (Aj^) = 0), can be written as 
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where a = sin "^[fcy/(e + 1)], a' = sin "^[fc.y/(e — 1)], fc = 

a is the angle of inci- 
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dence of the incoming electron (with wave- vector (fc, fcj,)) 
and a! is the angle of reflection of the Andreev reflected 
hole (with wave- vector (fc', fc^)). For retro AR, a', fc' have 
opposite signs from a, fc whereas for SAR, they have the 
same signs. The change from retro (e < 1) to SAR (e > 1) 
occurs at e = 1 (in our dimensionless units). 

Similarly for the superconducting barrier regions, the 
four component spinor solutions (u, v) contain electron 
wave-functions u of one valley and hole wave-functions 
V of the other valley. The DBDG equation can now 
be solved for any arbitrary energy e and the four solu- 
tions inside the superconducting barriers are given in the 
preprint version of Ref. |S| 
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where the subscripts 1/2 refers to the upper and lower 
signs on the RHS respectively, and similarly for 3/4 and 
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the barrier^. Moreover, for normal incidence, discrete 
Andreev bound levels are al so found between two super- 
conductors in graphene^^ESl. These levels can clearly lead 
to resonant transmissions in a SDB structure in graphene. 

We now follow the procedure set up in Ref. [25] and 
write a 4 X 4 S-matrix whose elements denote the var- 
ious net reflection and transmission amplitudes, except 
that now the electron and hole wave- functions, and con- 
sequently, the transmission and reflection matrices are 
each of them also 4-component. Hence, the §-matrix for 
the SDB structure in graphene for an incident electron 
with energy e can be written as 



§.= 



Mc Mac T, Tac 
T, Tac Kc Rac 

Tac Te ^Ac K, 



(10) 



where Kc stands for normal reflection of electrons or holes 
and M.Ac represents AR (and SAR), which is the reflection 
of an electron as a hole or vice-versa, from both the bar- 
riers combined. Similarly, Tc represents CT or normal 
transmission amplitude of electrons or holes while T ac 
corresponds to the non local CAR (and SCAR) amplitude 
for electron to hole conversion across the SDB structure. 
As mentioned earlier, each of these amplitudes are them- 
selves 4x4 matrices. The numerical results for these am- 
plitudes can be obtained by matching the electron and 
hole wave-functions between the normal and proximity 
induced superconducting graphene at each of the four 
interfaces and the numerical results for the resonances 
are discussed in Sec. IVl 
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Here, we have not taken the limit U ^ A,e. We have 
also obtained the solution for both right-moving and left- 
moving electrons and holes. 

Our aim now is to obtain the net quantum mechanical 
amplitudes for reflection, transmission, AR (and SAR) 
and CAR (and SCAR) of an electron incident on the 
SDB structure, after it has traversed both the barriers. 
A double barrier structure can always lead to resonances 
and this can affect the transmissions and the reflections 
through the system. For non-rel ativis tic electrons, this 
scenario has been studied befor^^iiasi p^j. relativistic 
electrons, the standard paradigm is that one cannot ob- 
tain conflned carrier states for normal incidenceSEl due 
to Klein tunneling. However, discrete energy levels can 
be found for carriers in graphene based quantum wires, 
as long as they have a non-zero component parallel to 



III. 



ANDREEV BOUND LEVELS 



Andreev bound states are formed due to multiple An- 
dreev reflections (for non zero incidence angle, by multi- 
ple retro andreev reflections only). Let us flrst consider 
the AR matrices which converts electrons to holes at each 
interface. The elements of the AR matrix can be derived 
using 
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Each of this conditions gives two equations, and we have 
only 8 non zero elements in the matrix M.a if there is no 
reflection. More explicitly, the matrix elements (M^)^ ., 
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Note that a needs to be small in order to ensure that the 
reflection is mainly retro-reflection. We will consider the 
effect of specular reflection (along with SAR) in the next 
section. In our analysis we have used left-right symmetry, 
but we have been careful to maintain the distinction be- 
tween electron and hole parameters, since exact electron- 
hole symmetry only exists at e = 0. At any finite energy, 
the symmetry is broken. We can also write down the 
4x4 transmission matrix T^ which allows for both nor- 
mal transmission and CAR through a single barrier. This 
matrix is defined by 



The other elements of the matrix M^ are found by in- 
terchanging a O —a', rA,he — >■ 'rA,eh and also changing 
their positions to the other off-diagonal block. In partic- 
ular, for a = 0, the matrix is much more simple and can 
be written as: 
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and we find that the matrix elements are given by 
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Clearly T^ includes both normal transmission and 
CAR. It is also clear that T^ can be written as T^ = 
T + T^ where T is a 4 x 4 matrix with two non-zero 
2x2 diagonal blocks and Ta is a 4 x 4 matrix with two 
non-zero 2x2 non-diagonal blocks. 

The phase matrix relating the electron and hole wave- 
function when it traverses the normal graphene region 
through a distance L (the distance between the two su- 
perconductors) is given by 
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a diagonal matrix with the entries 
'' ^) denoting the phases picked up 
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by the left and right moving electrons and holes respec- 
tively. The A and A' matrices which rotate the momen- 
tum operator to an arbitrary basis are given by 
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The condition for resonance or for a bound state in the 
normal graphene region between the two superconduc- 
tors is now just the condition that the total transmission 
computed as 



= TbM 
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has a vanishing denominator. This is precisely the con- 
dition for the Andreev bound states formed by retro ARs 
without the presence of reflection. From this condition, 
one can find the corresponding Andreev bound state en- 
ergy levels. Note that if we want the total transmission 
of electrons, we just need to replace T^ by T and if we 
want the total CAR of holes, we need to replace T^ by 
Ta- Both of them show the effect of the resonances. 



IV. EFFECT OF REFLECTION AND SAR 

For normal incidence in graphene, due to pseudo-spin 
symmetry, normal reflection is prohibited and we have 
pure Andreev bound states between the two supercon- 
ductors. Normal bound states (formed by multiple ordi- 
nary reflections) are, in any case, not possible in graphene 
even at any other incident angle, when ordinary reflection 
is allowed. This is because multiple ordinary reflections 




FIG. 2: (Color online) The electron and hole paths contribut- 
ing to the formation of Andreev bound levels between the 
two superconductors. The electrons have been shown as red 
lines, the retro AR, CAR holes as blue lines and the SAR, 
SCAR holes as green lines. The bound levels formed by mul- 
tiple retro AR have been shown as thick red and blue lines. 
However, not all possible paths have been shown in the figure. 
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FIG. 3: (Color online) The bottom graph shows the behaviour 
of angle resolved conductance, obtained numerically, in units 
of Ae^ /h as a function of energy in the subgapped regime (e ^ 
A) for ky = 0.75. Here, /:\/Ep = 0.05 and Uo/Ep = 10.0. 
The top graph depicts the denominator of Eq |20| for the same 
parameter values. 



are specular in nature and lead to a mode running along 
the y-axis, rather than a bound state. This is also true 
even if we have SAR, which leads to a specular Andreev 
mode running along the y-axig^l. 

However, at a non-zero angle of incidence, if we have 
specular reflection present (Andreev as well as normal) 
at each graphenc NS interface, we can still have Andreev 
bound states formed by multiple retro AR in between the 
two superconducting barriers. However, due to the spec- 
ular nature of the reflection and the SAR, they are no 
longer localised. This is shown in Fig. [2] The incident 



electron can transmit through the first superconducting 
barrier as an electron and then have multiple retro AR. It 
can also have either ordinary reflection or SAR from the 
interface. The transmitted electron once again (in fact, 
for several times, depending on the length of the graphene 
sheet in the j/-direction) can have multiple retro AR be- 
fore it is finally transmitted (or reflected) as an electron 
or a hole through the graphene SDB structure. The elec- 
tron can also cross the first superconductor as a hole by 
retro CAR and then have multiple retro-reflections from 
the two interfaces and then this process can continue as 
well. The number of reflections (or SAR 's) that can oc- 
cur in a given sample is controlled by the length of the 
graphene sheet in the j/-direction and the angle of inci- 
dence of the electron. Naively, the number should go as 
Ly/Lsina, where Ly is the length of the sheet in the 
y-direction and L and a have already been defined ear- 
lier. In Fig. [2] we have only shown some of the possible 
quantum mechanical paths to emphasize how the An- 
dreev bound states can form between the two barriers. 
We have also shown only a single incident electron, but 
the incident electron can also be at any point along the 
y-axis. Hence, if we measure the total output current 
collected throughout the j/-length of the graphene sam- 
ple, we should to able to get the signature of the many 
Andreev bound states present between the two supercon- 
ducting barriers. In fact, the signatures of the Andreev 
bound states should be present in all the four ampli- 
tudes Tc, rAc, tc and tAc, where these four amplitudes 
denote the quantum mechanical amplitudes for reflec- 
tion, AR (SAR), transmission and CAR (SCAR) across the 
SDB structure for an incident electron. They should show 
either a maximum or a minimum precisely at the Andreev 
levels obtained in Eq. [20j This is the main point that we 
wish to emphasize in this paper. The Andreev levels can 
be obtained by measuring the transmissions and reflec- 
tions through a SDB system even for large enough Ly-, 
it does not require an effective 'one-dimensional' system. 
Furthermore, the signature of the resonances occurs in 
all the four amplitudes. 



To illustrate the above argument, in FigjSJ we com- 
pare the positions of the resonances as a function of e/A 
obtained by the vanishing of the denominator of Eqj20j 
the analytically obtained Andreev levels, (which do not 
include the effects of reflection) with the numerically ob- 
tained values of the conductance (by solving the scatter- 
ing problem in the next section) for the same parameter 
values. Note that the numerical results do include reflec- 
tion since we have chosen ky = 0.75. Moreover, the nu- 
merical results include transmission throughout all points 
along the y-axis (width of the graphene sheet). But as 
can be seen from the figure, the resonances still fall on 
top of each other. This clearly shows that the presence 
of specular reflection (and SAR) has no effect on the po- 
sition of the resonances which simply occur due to the 
formation of the Andreev bound states discussed in the 
earlier section. 
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FIG. 4; (Color online) The behaviour of all possible quantum mechanical scattering probabilities (|rc|^, |tcl^, J^acI^, |tAcP) 
through the graphene SDB structure are plotted as a function of energy in the subgapped regime (e <C A) for three different 
values of ky. In (a), (b) and (c) solid red, green, blue and dashed black lines correspond to the |rAc|^, |'"cp, |iAcP and \tc\^ 
respectively. Here, A/Ef = 0.05 and Uq/Ef = 10.0. 



V. NUMERICAL RESULTS 

In this section we describe the consequences of all 
the allowed quantum mechanical processes across the 
SDB geometry in graphene. To the left of the SDB struc- 
ture, with an incident electron from the left, the wave- 
function can be written as 



V''=+ + Tc^- + rAcV 



(21) 



and to the right of the SDB structure, the wave-function 
can be written as 



icV'"^ + tAc'4> 



h+ 



(22) 



Here we use the standard wave-function matching tech- 
nique to solve such scattering problems to obtain all the 
four quantum mechanical amplitudes. Hence matching 
the wavefunctions for the normal and proximity induced 
superconducting regions (Eq.(4]l7)) at the four NS inter- 
faces in graphene {x = 0, a, a + L, 2a + L) forming the 
SDB structure, we obtain sixteen linear equations. Nu- 
merically solving these sixteen equations, we obtain the 
four amplitudes Tc, rAa tc and Iac, for the SDB struc- 
ture, for an incident electron with energy e below the 
gap A. Note that we distinguish between electron and 
hole parameters, and hence the four amplitudes will be 
different for incident electrons and holes. 

In our numerical analysis we do not distinguish be- 
tween the specular and retro Andreev reflections and we 
also allow for normal reflection at each NS interface in 
graphene, besides normal transmission and CAR (specu- 
lar and retro). The numerical results clearly show that 
for normal incidence of electron (a = 0) , the net normal 
transmission (CT) tc and the net AR rAc show resonant 
behaviour in the subgapped (e ^ A) regime. This is 
shown in the first panel in Fig. [4] and is precisely what 
one would expect, because for normal incidence Tc = 
due to the pseudo-spin symmetry, and this also leads 
to Iac = 0. For normally incident electrons, this kind 
of resonant behaviour is completely forbidden in a nor- 
mal double barrier (DB) structure in graphene due to 
the phenomenon of Klein tunneling. This is the striking 
difference between a normal DB and a SDB in graphene. 
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FIG. 5: (Color online) (a) The behaviour of the transmission 
resonances {T = \tc\^) is plotted as a function of energy in 
the subgapped regime (e <gi A) for three different values of 
a/L ratio. The red, blue and green lines correspond to the 
three different values of a/L which are 0.012, 0.017 and 0.022 
respectively, (b) The distance between consecutive resonances 
for the same three different values of a/L. For both the figures 
ky = 0.125. 



We also vary the momentum ky or equivalently the an- 
gle of incidence, and study the resonances in Fig. [3J As 
we have already mentioned, at fcy = or normal inci- 
dence, we have non-zero values only for tac and ic, and 
strong resonant behaviour for the transmission tc- As 
soon as the angle of incidence changes, we see the evo- 
lution of the tcTc^tAc and rAc in the panels in Fig. [3] 
As we increase the angle of incidence, due to finite fcy, 
normal reflection re between the barriers increases, and 
due to the presence of both re and rAc between the bar- 
riers the amplitude of the transmission resonances de- 
creases. However, at very large ky, the roles of reflection 
and Andreev reflection switch; this is true even for a sin- 
gle NS junction in graphene, since rAc is proportional to 
cos a and re is proportional to sin d^. Hence, for large 
ky, we find large values of r^, but very small values of rAc 
and once again, strong transmission resonances emerge. 

These resonances can also be tuned by varying param- 
eters such as the ratio of the width of the superconductor 
a to the length between the two superconducting barriers 
L. As a/L increases, we find that the resonances become 
sharper, and the distance between consecutive resonances 
increases. The behaviour of the resonances as a function 
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FIG. 6: (Color online) The behaviour of the angle-resolved differential conductance in units of Ae^ /h as a function of energy in 
the subgapped regime (e <^ A) for three different values of ky. Here, A/Ef = 0.05 and Uq/Ef = 10.0. 



of a/L is shown in Figs l4l^a) and WCb) . 

The net angle-resolved differential conductance 
through the SDB system is now given by 



G^N*Go \tAf cos(a') - \tf cos(a) 



(23) 



where N is the number of input channels or transverse 
modes in a graphene sheet of width Ly, Gq = 4e^//i is the 
unit of conductance and the factor of 4 comes from the 
pseudo-spin and valley degeneracies present in graphene. 
In our numerical analysis we have considered the tem- 
perature to be zero and also assumed linear bias. We 
again wish to emphasise the fact that since multiple An- 
dreev reflections between the two superconducting barri- 
ers not only include both specular and retro AR, but also 
normal reflections, the exit point from the second super- 
conductor can be anywhere along the width (j/-axis) of 
the graphene sheet. Therefore, both normal reflection as 
well as SAR can change the position of the transmitted 
beam along the y-axis at each reflection. Hence, the to- 
tal transmission here includes transmission at all points 
along Ly. In other words, the output lead has also to 
be as wide as the graphene sheet. In Fig. [5} we show 
the net angle resolved Landauer-Buttiker conductance, 
given in Eq. [231 (for A^ = 1) as a function of the energy 
of the incident electron e ^ A, where again, the reso- 
nant behaviour can be seen for different values of ky or 
the incident angle. The behaviour of the conductance, 
also shows how the Andreev levels evolve as a function of 
the incident angle, showing that the height of the reso- 
nances is large when the multiple reflection between the 
barriers is either dominated by retro AR (small angles) 
or normal reflection (large angles). 



VI. 



SUMMARY AND DISCUSSIONS 



To summarize, in this paper, we have computed the 
transmission of an electron through a SDB structure in 
graphene and shown the resonant suppression of Andreev 
reflection at certain energies below the superconducting 
gap A where normal transmission \tc\'^ becomes unity. 
This resonant behaviour is absent in a normal double 
barrier in graphene due to Klein tunneling. We also show 



that the resonant suppression is due to the formation of 
Andreev bound states between the two superconducting 
barriers. Even at finite incident angles, the position of 
these Andreev bound levels remain unchanged in pres- 
ence of reflection and SAR. Although, the transmission 
resonances get damped as the incident angle increases 
due to reflection and SAR. However, at large angles, the 
roles of AR and ordinary reflection get reversed, and once 
again, we see strong transmission resonances. The point 
that we wish to emphasize is that at any angle of in- 
cidence, all four transmission and reflection amplitudes, 
as well as the total conductance, shows strong signals of 
the Andreev bound states. We have also studied the res- 
onances as a function of the ratios of the width of the 
superconducting barriers and the length of the normal 
graphene region between the two superconducting barri- 
ers. The most interesting study is the evolution of the 
resonances as a function of ky, the momentum in the 
y-direction or equivalently the angle of incidence of the 
incoming electron. 

As far as the practical realization of such a SDB struc- 
ture in graphene is concerned, it should be possible to 
fabricate such a geometry by depositing thin strips of a 
spin singlet superconductor (like Al or Nb) on top of a 
graphene sheelP at two places. The width of the strips 
should be of the order of the superconducting phase co- 
herence length (10 — 15nm in case of Nb) for CT and 
CAR (and SCAR) to take place. For a given a/L ratio, the 
|icP — 1 resonance in this SDB geometry can be tuned 
by varying the energy of the incident electron (which can 
be done by applying a small bias voltage between the 
two reservoirs keeping within linear response, so that our 
calculations are valid). In Fig. l4[a) and 6(a) the equiv- 
alent temperature at which the flrst resonance occurs is 
approximately lOmK (critical temperature Tc ~ 9.2K in 
case of Nb) for a ~ 10 — 15nm and L ~ 2^m. 

Note also that, experimentally, it is not easy to con- 
trol the angle of incidence of the electrons that impinge 
on the SDB structure. In fact, for a given sample what 
can possibly be done is to grow the barriers at different 
tilted positions by lithographic techniques which allows 
the electrons to be incident on the barriers at different 
angles. However, in a given experiment, the angle of 
the incident electron is flxed, since the incident electrons 



are collimated. Therefore different experiments need to 
be performed with baUistic graphene samples to obtain 
results at different angles of incidence. In a given experi- 
ment, the angular spread of the incident electrons which 
would exist due to lack of perfect coUimation is actu- 
ally quite small, hence even after angular averaging our 
results at a given angle should qualitatively remain the 
same. 

We also expect the features of the resonances to be 
qualititavely unaffected by long-ranged impurities (slow 
on the scale of the lattice spacing i.e. kpl^ ^ 1 where 
kp ^ l/d, d is the lattice spacing and Z„i is the mean- free 
path of electron in graphene), because such impurities 
in graphene can only cause intra-valley scattering. On 
the other hand, short range impurities {kpl^ <§C 1) can 
cause intervalley scattering which can destroy the res- 
onance. Also, strong disorder can localize electron/hole 
states between the two barriers which always can destroy 
the resonance. Another experimental variable which is 
expected to destroy the resonance is the presence of a 
magnetic field which bends the electron paths between 
the two barriers. Hence, it prevents the multiple retrac- 
ing of the path, which is needed for resonance. Therefore, 
we do not expect resonances in the presence of strong dis- 
order or even weak magnetic fields. 

In this work, we have restricted ourselves to rectan- 
gular barriers, whereas experimentally, the barriers are 
expected to be smootheil^Sl. As shown in the references 
by Cayssol et al.^'^ and SonirpS, the transmissions are af- 
fected by the shape of the potential, and it is shown that 
the resonances for normal double barriers in graphene are 
sharpened by smoother barriers. By analogy, we expect 



our resonances to be sharpened if we make the barriers 
smooth instead of using rectangular barriers. Neverthe- 
less, to check it explicitly, we need to repeat the calcu- 
lations for other shapes of barriers which is beyond the 
scope of this paper. 

From the application point of view, this SDB geometry 
also can be used to confine electrons in graphene which 
is experimentally very challenging due to Klein tunnel- 
ing. The resonant structure of |icP as a function of e/A 
is analogous to the Coulomb blockade peaks of a quan- 
tum dot. Hence, this kind of SDB geometry can be used 
to form quantum dots in graphene. Another interest- 
ing point which can be mentioned is that this kind of 
proximity induced SDB geometry can also be obtained 
by replacing the graphene sheet by a 2D Quantum Spin 
Hall topological insulator. The only difference will be the 
replacement of the factor of 4 in the Landauer-Buttiker 
conductance formula by unity, since there are no pseudo- 
spin or valley degeneracies present in such quantum spin 
Hall insulator. 



Acknowledgments 

The work of A.S. was supported by the Feinberg fel- 
lowship programme at WIS, Israel. One of us (A.S.) ac- 
knowledge Anindya Das for useful discussions and com- 
ments. One of us (S.R.) would also like to acknowl- 
edge hospitality at the Perimeter Institute for Theoretical 
physics, Canada, where this work was completed. 



^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. 

Firsov, Science 306, 666 (2004). 
^ A. K. Geim and K. S. Novoselov, Nat. Materials 6, 183 

(2007). 
^ A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 

Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 

(2009). 
* S. D. Sarma, S. Adam, E. H. Hwang, and E. Rossi, 

Electronic transport in two dimensional graphene (2010), 

arXiv: 1003.4731 [cond-mat.mes-hall]. 
^ M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, Nat. 

Phys. 2, 620 (2006). 
^ H. B. Heersche, P. J. Herrero, J. B. Oostinga, L. M. K. 

Vandersypen, and A. F. Morpurgo, Nature 446, 56 (2007). 
'' A. F. Andreev, Sov. Phys. JETP 19, 1228 (1964). 
® C. W. J. Beenakker, Phys. Rev. Lett. 97, 067007 (2006). 
® C. W. J. Beenakker, Rev. Mod. Phys. 80, 1337 (2008). 
^° G. E. Blonder, M. Tinkham, and T. M. Klapwijk, Phys. 

Rev. B 25, 4515 (1982). 
^^ S. Bhattacharjee and K. Sengupta, Phys. Rev. Lett. 97, 

217001 (2006). 
^^ S. Bhattacharjee, M. Maiti, and K. Sengupta, Phys. Rev. 

B 76, 184514 (2007). 
^^ J. Linder and A. Sudb0, Phys. Rev. Lett. 99, 147001 



(2007). 
" J. Linder and A. Sudb0, Phys. Rev. B 77, 064507 (2008). 
^^ M. Titov and C. W. J. Beenakker, Phys. Rev. B 74, 

041401 (R) (2006). 
^^ M. Maiti and K. Sengupta, Phys. Rev. B 76, 054513 

(2007). 
^^ D. Greenbaum, S. Das, G. Schwiete, and P. Silvestrov, 

Phys. Rev. B 75, 195437 (2007). 
^* J. Cayssol, Phys. Rev. Lett. 100, 147001 (2008). 
'^^ C. Benjamin and J. K. Pachos, Phys. Rev. B 78, 235403 

(2008). 
2" V. V. Cheianov and V. I. Falko, Phys. Rev. B 74, 

041403(R) (2006). 
^^ J. M. Pereira, P. Vasilopoulos, and F. M. Peeters, Appl. 

Phys. Lett. 90, 132122 (2007). 
^2 C. Bai, Y. Yang, and X. Zhang, Physica E 42, 1431 (2009). 
^^ C. W. J. Beenakker, A. R. Akhmerov, P. Recher, and 

J. Tworzydio, Phys. Rev. B 77, 075409 (2008). 
^^ A. F. Morpurgo and F. Beltram, Phys. Rev. B 50, 1325 

(1994). 
'^^ A. Kundu, S. Rao, and A. Saha, Euro. Phys. Lett. 88, 

57003 (2009). 
'^^ A. F. Young and P. Kim, Nat. Phys. 5, 222 (2009). 
^^ J. Cayssol, B. Huard, and D. Goldhaber-Gordon, Phys. 

Rev. B 79, 075428 (2009). 



E. B. Sonin, Phys. Rev. B 79, 195438 (2009). 



